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Abstract 

A higher grading afflne algebraic construction of integrable hierarchies, contain- 
ing the Wadati-Konno-Ichikawa (WKI) hierarchy as a particular case, is proposed. 
We show that a two-component generalization of the Schafer- Wayne short pulse 
equation arises quite naturally from the first negative flow of the WKI hierarchy. 
Some novel integrable nonautonomous models are also proposed. The conserved 
charges, both local and nonlocal, are obtained from the Riccati form of the spectral 
problem. The loop-soliton solutions of the WKI hierarchy are systematically con- 
structed through gauge followed by reciprocal Backlund transformation, establishing 
the precise connection between the whole WKI and AKNS hierarchies. The connec- 
tion between the short pulse equation with the sine-Gordon model is extended to 
a correspondence between the two-component short pulse equation and the Lund- 
Regge model. 

1 Introduction 

The classification of integrable models is of fundamental importance. A systematization in 
this vein is the graded afhne algebraic construction of integrable hierarchies |HI2], unifying 
in a single algebraic structure, different models that at first sight appear unrelated. For 
instance, the modified Korteweg-de Vries (mKdV) and sine/sinh-Gordon models are mem- 
bers of the same hierarchy. Another example is given by the cubic nonlinear Schrodinger 
(NLS) and the Lund-Regge models, which are members of the Ablowitz-Kaup-Newell- 
Segur (AKNS) hierarchy. Due to the fact that all the models within a hierarchy share 
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the same algebraic structure, solitonic solutions for all of them can be constructed in a 
universal manner, through the dressing method [3J. 
The short pulse equation (SPE) 



\ u xt — u 



+ 1 (« 3 L (i-i) 



was proposed [4J as a model to describe ultra-short optical pulses traversing within a non- 
linear media. This equation replaces the standard NLS equation in the ultra-short pulse 



regime, being a very good approximation of Maxwell's equations. Since then, (1.1) was 



proved to be integrable [5j having a Lax pair of the Wadati-Konno-Ichikawa (WKI) type 



[6]. Moreover, it was shown [5] that (1.1) is related, through a hodograph transformation, 
to the sine-Gordon equation 

e ys = 4sm8 (1.2) 

from whence its solitary solutions were first obtained [7]. The hodograph transformation 
was further explored [8J to construct solitonic and quasi-periodic solutions. The transfor- 



mation connecting (1.1) and (1.2) reads 



dx = ady — 2u 2 dt, 
dt = ds. 



-i (1-3) 
a 2 = (l + u 2 ,) = cos 9, 

u = \e s . 



Through a recursion operator approach, it was proposed in [9] a hierarchy containing ( 1.1 ) 
and also the elastic beam equation (EBE) [10J 
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The bi-Hamiltonian structure and nonlocal conservation charges for this single component 
field hierarchy were also proposed [llj. Recently, a multi-component generalization of the 
short pulse equation was proposed [12j together with its soliton solutions. 

On the other hand, the algebraic construction [Tj |2] of integrable hierarchies is based 
on Toda field theories and the fields appear as coefficients of a linear combination of zero 
grade operators. This construction was further generalized \Vi\ IT4"1 ITS] by the addition of 
fields associated to higher grading operators, yielding the generalized affine Toda models. 
The higher grade fields are physically interpreted as matter fields with the usual Toda 
fields coupled to them. 

In this paper, we propose a general higher grading construction for the zero curvature 
equation, containing the WKI hierarchy as a particular case. In our construction, the 
zero grade Toda fields are completely removed, remaining the higher grade fields only. 
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Considering the WKI hierarchy, its previous known models appear as members of the 
positive (time) flows. We then extend the WKI hierarchy to incorporate negative flows, 
which were not been previously considered. Deriving the first negative flow, a model which 
is a two-component generalization of the SPE (1.1) emerges quite naturally, showing that 
the SPE is in fact a natural member of the WKI hierarchy and underlies the same algebraic 
structure as the previous known models within the positive flows. This two-component 
short pulse equation (2-SPE) was also introduced recently in [12J under a different context. 
We also introduce some novel nonautonomous models, mixing a positive flow with a 
negative one [IE]. The conserved charges, both local and nonlocal, are constructed from 
the Riccati form of the spectral problem [T7] . 

The WKI and AKNS hierarchies are known to be related [T8| IT9] . but this connec- 
tion was not sufficiently explored. Driven by this results we establish, from an algebraic 
perspective, the relation between both hierarchies, thus systematizing and extending the 
hodograph transformation (1.3) for the whole WKI hierarchy. As an example, we explic- 
itly derive a correspondence between the 2-SPE and the Lund-Regge model, where (1.3) 
appear as a particular case. Our relations also provide explicit and systematic solutions 
for all the models within the WKI hierarchy in terms of the known AKNS tau functions. 

We organize our work as follows. In section [2] we introduce the main algebraic con- 
cepts and the usual construction. In section [3j a higher grading integrable hierarchy with 
an arbitrary affine Lie algebra Q is proposed, then, in section [ij this construction yields 
the WKI hierarchy by choosing Q = A\ with homogeneous gradation. The known mod- 
els of the WKI hierarchy are derived by considering the positive flows and the 2-SPE 
emerges from the first negative flow. We also propose some novel nonautonomous inte- 
grable models. In section |5j we derive the local and nonlocal conserved charges for the 
WKI hierarchy. In section |6j we show that the dressing method is unable to solve the 
WKI hierarchy, then we gauge the AKNS hierarchy into the WKI hierarchy with the aid 
of a reciprocal Backhand transformation, mapping the corresponding flows of both hier- 
archies. We consider explicitly this mapping between the first negative flows, obtaining 
a correspondence between the 2-SPE and the Lund-Regge model. Our final remarks are 
presented in section [7} 



2 The usual algebraic construction 

Let Q be an affine Kac-Moody algebra and Q an operator decomposing the algebra into 
the graded subspaces 

g = Q)gv\ [Q,g {j) ]=3Q {3) . (2.1) 



3 



As a consequence of the Jacobi identity, [g {i \G {j) ] C . Let E be a semi-simple 

element, with a definite grade, defining the kernel subspace 

K = {T G G I [S,T] = 0}. (2.2) 

The image subspace, At, is its complement and Q = /C © At. Then we have the relations 
\JC, fC] C /C, [/C, At] C At and we assume the symmetric space structure [At, At] C /C. 
Integrable hierarchies are constructed from the zero curvature equation 

[d x + U,d t + V] =0, (2.3) 

where the Lax pair U and V lie on Q and have the following form, 

U = + A (o) [0] (2.4) 

where E = E^ G is a constant semi-simple element, defining /C and At, and A^ o ^[0] G 
At*- -*, where At*- - 1 — At H f^ -*, is the operator containing the fields of the models. For 
instance, if At (0) is spanned by {jj®, . . . , T^ 0) } then A( o )[0] = 1 T 1 (O) + ••• + 0^ o) , 
0j = 0j (x, t) and <p = (</>!, . . . , 0h) T - The other Lax operator has the form 

n 

V = ^2D {i) (for positive flows), (2.5) 

i=0 

-1 

y = ^ ,D {i) (for negative flows). (2.6) 



The positive integer n labels the (time) flow under quest. Formally, we should have written 
V = V (n) and t = t n , showing explicitly the flow we are considering, i.e. the particular 
model within the hierarchy, but we usually omit the index n that will be clearly indicated 
from the context. It is possible to show that the operators G are determined 



in terms of the fields contained in by projecting the zero curvature equation (2.3) 



into each graded subspace according to (2.1). This procedure enables us to construct the 



operator V, one for each flow n, in terms of the algebraic structure of the operator U, thus 
yielding a nonlinear model together with its Lax pair. Therefore, by construction, all the 
models are immediately classified through (Q,Q, E^), defining an integrable hierarchy of 
differential equations. 



3 A higher grading construction 

We now propose a higher grading construction without the usual Toda fields, keeping the 
higher grade fields only. Let us define the following zero curvature Lax pair 

U = E^ +A^[(f>], (3.1a) 
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n 



v=Yl D(i) M> ( 3 - lb ) 

i=—m 

where E = G is a constant semi-simple element and G .M 1 ' 1 *', where 

= M. R is the operator containing the fields now having grade one. The 

positive integers n and m label each model, mixing a positive flow n with a negative flow 
m |16j . To study positive and negative flows separately one should consider 

V = ^2D {i) (for positive flows), (3.2) 

i=l 
1 

V = D(i) ( for negative flows). (3.3) 

i=—n 

With this algebraic structure, the zero curvature equation can be solved non trivially in 



the following way. Taking (2.3) with (3.1 ), the projection into each graded subspace yields 



the following set of equations 

[ E W + £>(»)] = grade n + 1, 
d x DW + + A®, D^V] = grade n, 

d x DV> - d t A^ + [£« + D<®] = grade 1, (3.4) 

d x D { ~ m+1) + + A {1 \ D { ~ rn) ] = grade -m + 1, 

d x D ( - m) = grade -m. 

Note that the equations of motion comes from the projection into AidQ^. Each equation 
is further decomposed into JC and A4 components. For mixed or negative flows, we have 
the equations of motion in the following form 

dtAV = d x D% + [E<»,D$l + [A« £>£>], (3.5) 

while for positive flows the equations of motion are given by 

cU (1) = d x D$, (3.6) 

where = D^' +D^. We remark that this construction is valid for an arbitrary affme 
Lie algebra Q with a grading operator Q and arbitrary grade one semi-simple element 
EW. 
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4 The Wadati-Konno-Ichikawa hierarchy 

Let us now consider the Kac-Moody algebra A\ = {if- 7 ,E+, E 3 _ ,c,d}, where j is an 
integer, with commutation relations given by 

[H\ H 3 ] = 2k5 k+j , c, [H\ E{] = ±2El +j , [E k + , El] = H k+3 + kS k+j , c, 
[d,T j ] = jT 3 , [c,T 3 ] = 0, where T j e {H j , E{, El) . 

In the construction of integrable models we use the loop algebra only, achieved by setting 
c = 0. The homogeneous gradation Q = d yields the grading subspaces 



?00 



{Hi,E{,Ei}. (4.2) 



Let us fix the semi-simple element as E^ = H 1 , yielding KS 3 ^ = {H 3 } and M.^' = 
{E 3 + , E 3 S\- Therefore, the operator containing the fields must have the form 

= q(x, t)E\ + r(x, t)E x _, (4.3) 



where q and r are both fields of the models. The Lax operator V in (3.1b) is a sum of 



elements in the form D^' = a,jE 3 + + bjE J _ + CjH J , whose coefficients will be determined 



in terms of q and r by solving the zero curvature equation as described in (3.4). 
4.1 Positive flows 



The positive flows (3.2) are constructed from the zero curvature equation 

[d x + H 1 + qE\ + rEl, d t + D {n) + + • • • + L> (1) ] = 0. (4.4) 
Now let us consider the first models. 

n = 1. Yields the trivial equations q t = q x and r t = r x . 

n = 2. Starting from the projection we have a 2 = c 2 q and b 2 = c 2 r. The Q^ 1 
projection gives 

ai = ciq - \d x a 2 , h = c x r + \dj> 2 , d x c 2 = ra x - qbi, (4.5) 
from which we have 

d x lnc 2 = —\d x In (1 + qr) , c 2 = (1 + qr)~ l/2 . (4.6) 
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The K, component of the projection implies d x C\ = 0. Choosing c\ = c[j] the Ai 
component yields the following equations of motion 



d t q + d 2 x 



2(1 + qr) 



1/2 



0, d t r - dl 



2(1 + qr) 



1/2 



0. 



(4.7) 



From the coefficients a,, b$ and q just found we have the explicit Lax pair for this model 
which reads 



U = H 1 + qE]+ rE l _ 



V 



Q 



:E 2 _ + 



[1 + qr) 1 ^ 2 + (1 + qr) 1 ^ 2 (1 + qr) 1 ^ 

q \ _i „ / r 



H 



2(1 + qr) 



1/2 



Ei + d x 



2(1 + gr) 



1/2 



(4.8a) 



(4.8b) 



The choice r = q*, x — >■ zx and t — >■ 2it reduces (4.7) to the Wadati-Konno-Ichikawa- 
Shimizu equation |2Uj 
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idtq + di 



2^1/2 



0. 



(4.9) 



n = 3. In this case we get the following model 



$9 - <9 2 



qi 



4(l + qr) 3/2 
whose Lax pair is 

U = H 1 + qE\ + rE\ 
V 



0, d t r - d x 



4 (1 + qr) 



3/2 



(4.10) 



(4.11a) 



(1 + qr) 



" El + 



1/2- 



1 + qr) 



1/2- 



9i t^2 



Ei + 



2(l + gr) 3/2 + 2(l + gr) 3/z 4(l + gr) 



3/2 



1 + gr) 
E 2 _ + 



1/2- 

q x r-qr x r2 



,3/2 



+ 3* 



4(1 + gr) 



3/2 



sl + a 



4 (1 + qr) 



3/2 



s 1 . 



(4.11b) 



Choosing r = q and t — > —At the model (4.10) becomes 



dtq + ^ 



(1 + g 2 ) 



2^/2 



(4.12) 



1 Without loss of generality because it contributes with a linear term that can be absorbed by a simple 
change of coordinates. 
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which is the EBE (1.4). Note also that, by choosing r = 1, q — v — 1 and t — > —2t, the 
system (4.10) is transformed into the Dym equation 



(4.13) 



The Lax pair for the Dym equation, and also for (4.12), are obtained from (4.11) as 
particular cases. The models (4.7) and (4.10) were proposed in [6J and are the first 



members of the positive WKI hierarchy. 
4.2 Negative flows 

We now extend the WKI hierarchy to incorporate negative flows. According to our con- 



struction, the negative flows are generated using the operator (3.3) yielding 



[d x + H 1 + qE\ + rEl, d t + D^ n) + D^ n+1) + ■■■ + D w + D w ] = 0. 



(4.14) 



n = 1. The projection into ^ implies that a_i, 6_i and c_i are all constants^} The 
g(o) projection yields 

d x ao — 2c_ig — 2a_i, d x bo — — 1c-\T + 26_i, d x co = a^\r — b-iq. (4.15) 
Let us introduce the operator 

d^fix) = f X f(y)dy. (4.16) 



We assume that the fields and its derivatives of any order decay sufficiently fast when 
\x\ — > oo. Under this condition d x d~ 1 f{x) = d~ 1 d x f(x) = f(x). 



To eliminate an explicit x dependence in (4.15) we set a_i = 6_i = thus 
a = 2c_i<9~ 1 g, b = — 2c_i9~ 1 r, Co = const. 



(4.17) 



The projection into implies that a x = c^q and b\ = C\r. The projection yields 
the field equations plus one constraint 

d t q = d x ai + 2a - 2c q, d t r = d x b x - 2b + 2c r, d x c x = a r - b q. (4.18) 



From this last equation and (4.17) we have 



ci = 2c^d x 1 (rd x l q + qd x V) 



(4.19) 



2 In fact they can depend on t, but this contributes with an overall factor that can be absorbed by a 
change of coordinates, so we assume they are constants. This t dependency will be exploited to construct 
nonautonomous models in the following section. 
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Choosing Co = (j^jand fixing c_i = 1, we obtain the nonlocal equations of motion 



d t q = 4d x 1 q + 2d x (qd x 1 (rd x l q + qd x V)) , 
d t r = Ad~ x r + 2d x {rd~ l (rd^q + qd~ x r)) . 



Introducing the new fields defined by 

q = d x u, 



r = d x v, 



we can write (4.20) in a local form yielding 



u xt = 4m + 2d x (uvu x ) 
v xt = Av + 2d x (uvv x ) 



(4.20a) 
(4.20b) 



(4.21) 



(4.22a) 
(4.22b) 



This model was also proposed recently in |T2"|^| Note that (4.22) is symmetric by the 
interchange of u and v and invariant under the following transformations 



x — >■ /J.X, 



t — >• n H, 



ii — y n u, 



v ->■ [i 2 n v, 



for a constant (jl. The explicit Lax pair of model (4.22) is given by 



U = H 1 + u x E 1 + + v x E 1 _, 

V = H~ l + 2uE°, - 2vE°_ + 2UVH 1 + 2uvu x E\ + 2uvv x E 1 _. 



(4.23) 



(4.24a) 
(4.24b) 



If we reduce the model (4.22) to a single component field by choosing v — u, we have the 



SPE (1.1) 



\U xt 



u 



(4.25) 



Therefore, the model (4.22) is a two-component generalization of the SPE and a natural 



member of the WKI hierarchy, sharing the same algebraic structure as the previous known 
models within the positive flows. A n-component generalization of the SPE, following the 



structure of (4.22), was also proposed in |12j and our construction strongly suggests that 
it can be obtained by considering the untwisted algebra A n _\ ~ s£(n), which will also 



generate n-component extensions of (4.9) and (4.12). 



3 It contributes with a linear term that can be absorbed by change of coordinates. 

4 Choosing u = u + av and v = u — av where a — ^/=pl, after summing and subtracting the equations 
we get \u x t — u + ^d x ((w 2 ± u 2 ) u x ^j and \v x t = v + ~d x {{u 2 ± -0 2 ) v x ^. Under a small amplitude limit 
v <C 1 we have \u xt — u + g (u 3 ) and \v x t — v + \d x (u 2 v x ), which were recently considered in |21| 
and represents the small field v interacting with the background field u, being a free solution of the SPE. 
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n = 2. Setting c_2 = 6-2 = to eliminate explicit x dependencies and introducing new 
fields defined through q = u xx and r = v xx , we get the following model 



vu, 



))} 



u xxt = 2a {2u x + d x (u x v x u xx )} 

— 46 {2-u + u x v x u xx + d x [u xx (uv 

v xxt = 2a {2v x + d x (u x v x v xx )} 

+ 46 {2v + u x v x v xx - d x (v xx (uv x - vu x ))} 



where a = c_x and b = c_ 2 are constants. The Lax pair for this model is 



(4.26a) 
(4.26b) 

(4.27a) 



U = H 1 + u xx E\ + v xx E 1 _, 

V = bH~ 2 + 2bu x E- 1 - 2bv x EZ 1 + aH' 1 

+ 2 (au x - 2bu) El -2 (av x + 2bv) E°_ + 2bu x v x H° 

+ 2u xx (au x v x — 2b (uv x — vu x )) E + + 2v xx (au x v x — 2b (uv x — vu x )) E_ 

+ 2(au x v x -2b(uv x -vu x ))H 1 . (4.27b) 



Note that the term next to the a coefficient in (4.26) is the system (4.22). Setting a = 



and b = | we have 

(u xt + u xx (uv x - vu x )) x = -2u - u x v x u xx , 
(v xt + v xx (uv x - vu x )) x = 2v + u x v x v xx , 

which after the substitution v = u* and t — >■ it becomes the single equation 

(iu xt + u xx (u*u x - uu* x )) x = 2u + u x u* x u xx . 



(4.28a) 
(4.28b) 



(4.29) 



4.3 Small amplitude limit 

Let us give a physical motivation by considering a small amplitude limit. Consider the 



reduced second positive flow (4.9). Assuming that q is small we have 

(i+ki 2 r /2 =i-^i 2 +--- 



and then (4.9) becomes 



Consider a standard multiple scale expansion (e<l) 



q = ^2e n q n (x ,t ,x 1 ,t 1 , . . . 



U = e n t. 



(4.30) 



(4.31) 



(4.32) 



n=l 
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Taking order of e in (4.31) we obtain the plain wave solution 

Ql = Ae ia + A*e~ ia , a = kx - k\, A = A(x 1 ,t 1 , ...). (4.33) 

Taking order of e 2 , the elimination of the secularity requires that A = A(xi—2Kt±, X2, £2, • • • ) 
and we choose q2 = 0. Finally, the lower nonlinear contribution comes from the order of 
e 3 and after changing to a coordinate system moving with the envelope A, £j = X{ — 2nti 
and Ti = ti, we obtain the NLS equation 

id n A + dlA+— \A\ 2 A = 0. (4.34) 
Up to this order the small amplitude solution reads 

q = e {Ae iu + A*e~ ia ) + e 3 (Be 3la + B* e - 3i(T ) + O (e 4 ) , (4.35) 



where A is determined from (4.34) and contains a nonlinear contribution, and B comes 
from higher order terms. The positive flows of the WKI hierarchy contain a strong nonlin- 
earity and describes large amplitude solutions, which in the small amplitude limit recover 
the behaviour described by the AKNS models, like the NLS equation. Due to the large 
amplitude solutions these models may be of particular interest in nonlinear optics, plas- 
mas, Bose-Einstein condensates and water waves. 



Performing the same procedure for the EBE (4.12), the order of e gives the linear wave 



(4.33) with a = kxq + /t 3 to- The order of e 2 implies A = A {x\ + 3/t 2 ti, x%, ti, ■ ■ •) and 
q2 = 0. The order of e 3 , after changing to a coordinate system & = Xi + 3/t 2 tj and r, = ij, 
gives again the NLS equation 

id T2 A-3KdlA-3K 3 \A\ 2 A = 0. (4.36) 



Consider the SPE (4.25). The e order gives a linear wave U\ having the same form as 
(4.33) with a = kxq — 4/t _1 to- The e 2 order implies A = A (x\ + 4/t _2 ti, X2, ^2, • • •) and 
U2 = 0. The e 3 order, in the coordinates £j = x^ + An~ 2 ti and r< = tj, gives 

id T2 A + 4K~ 3 dlA + 2k\A\ 2 A = 0. (4.37) 

Note that the sine-Gordon equation u x t = smu ~ u — ^u 3 + - ■ ■ has stronger nonlinearities 
than the short pulse equation, but up to order e 3 both of them yields qualitatively the 
same multiple scale approximation. Contrary to the positive flows of the WKI hierarchy, 
the short pulse equation does not seem do describe large amplitude solutions. 

4.4 Mixed flows 

It is possible to combine a positive flow n with a negative flow m by considering the 
operator ( 3.1b[ ). Let n = 2 and m — 1, leading to the zero curvature equation 

[d x + H 1 + qE\ + rEl, d t + £> (2) + D w + D<® + D^] = 0. (4.38) 
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We solve each grade projection starting from highest to lowest, exactly as before, the only 
modification is that we can let some coefficients, that were previously considered constants, 
now to depend on t, thus providing the nonautonomous ingredient. For individual flows 
these coefficients were not interesting because they come as a global factor in the final 
equation. Therefore, after solving (4.38) we obtain 

'■9. 



Vxt 



a{t) a2 



u. 



(1 + u x v x ) 1/2 



4b(t) ( u + -d x (uvu x ) 



V.r 



+ 4b(t) ( v + -d x (uvv x ) 



(4.39a) 
(4.39b) 



(4.40a) 
(4.40b) 



(1 + u x v x ) 1 ' 2 

where a(t) and b(t) are arbitrary functions. Its Lax pair given by 

U = H 1 + u x El + v x E l _, 
V = a(t)V(+2) + b(t)V(-l), 

where V(+2) stands for the operator corresponding to the second positive flow ( |4.8b ), 
with q = u x and r = v x , and V(— 1) is (4.24b). Model (4.39) is a nonautonomous mixture 
of (4.7) and (4.22). This model also admits the reduction v = u* , t —> it and x — > ix, 

] -d 2 



yielding 



u x t 



ia{t) 2 



u x 



1 



12^/2 



-4b(t)[u--d x (\u\ 2 u x ) 



In the same way, for n = 3 and m — 1, we have 



Uxt 



Vxt 



fl ( f ) Q2 



Og) 2 

4 x 



u x: , 



XX 



[I + U X V. 



\3/2 



+ 



Ab{t) (u + -d x (uvu x ) 



x v x) 



[I + U x V x j 



\3/2 



+ Ab{t) (v + l -d x {uvvx) 



together with its Lax pair 



U = H l + u x E 1 + + v x E l _ ) 
V = a(t)V(+3) + b(t)V(-l). 



Choosing v = u, (|4.42|) leads to a combination of EBE and SPE 

Uxx 



a 



xt 



a(t) 2 
4 x 



^2^ 3 / 2 



+ 4b(t) [u + - (u 3 ) i 



(4.41) 

(4.42a) 
(4.42b) 



(4.43a) 
(4.43b) 



(4.44) 



;i + {ux) 

It is possible to combine any positive flow with a negative one, generating nonau- 
tonomous models inside the hierarchy. Due to a(t) and b(t) the dispersion relation will 
have a time depend velocity and the solitons will accelerate. This models may be nice 
candidates in applications having accelerated ultra-short optical pulses 



12 



5 Conservation laws 



Consider the linear spectral problem associated to the WKI hierarchy 

(d x + U)V = 0, (d t + V)V = 0, 



(5.1) 



where U = H 1 + qE+ + rE\ and V depends on the particular flow. Using a matrix 
representation we have 



U 



X Xq 
Xr -X 



V 



Vu(X) V 12 (X) 
V 2 i(X) V 22 (X) 



'02 



(5.2) 



where A is the loop algebra spectral parameter. From ( |5.1 ) we write the Riccati form of 
the spectral problem [T7], whose compatibility yields the conservation laws 



d t (qV) - d x (A- 1 (V n + V 12 T)) = 0, 
d t (rT- 1 ) - d x (A" 1 (V 22 + V^T- 1 )) = 0, 



where we have defined 



r = 



^2 
1>1 



(5.3) 
(5.4) 

(5.5) 



Therefore, we can use qT and rT^ 1 to construct an infinite number of conserved charges 
assuming a power series in A. Deriving T and with respect to x and writing the result 
in terms of F = qT and G = rT" 1 , we obtain the Riccati equations 

qd x (jj + Xqr - 2AF - AF 2 = 0, (5.6) 

rd x (^A + Xqr + 2AG - XG 2 = 0, (5.7) 

which are the generating equations for the conserved densities. Both equations lead to 
the same results, the densities obtained from (5.7) are the same as those from (5.6) with 
the interchange q •<->■ r, so we will focus on (5.6) only. 



5.1 Local charges 

Let us consider the expansion 

oo 

F = Y,fnX~ n . (5.8) 

n=0 
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Substituting (5.8) into (5.6) and taking powers of A we determine each f n recursively, the 



first ones are given by 



/o 
fi 



[1 + qr) 



1/2 



-d x In f — . 

2 I (l + gr) 1/2 



+ 



2q (1 + qr) 



1/2 ' 



(5.9) 
(5.10) 



The charges associated to these densities, which are conserved thanks to (5.3), are 

/oo 
(l + qr) 1/2 dx, 
-oo 



oo 



2q (1 + qr) 



1/2 



(5.11) 
(5.12) 



These charges coincide with those in [6] and are the Hamiltonians generating the models 
within the positive flows. 

5.2 Nonlocal charges 



The most general expansion compatible with the Riccati equation (5.6), able to generate 
nonlocal densities, is given by 



F = E f-^ n - 



(5.13) 



n=-l 



Substituting (5.13) into (5.6) and taking the lowest order, A 1 , we obtain the Bernoulli 

d x h-{d x \nq)h = fl 



differential equation 
whose solution is 



(5.14) 



f 1 = -d x \n(d- 1 q). (5.15) 

Note that it is a total derivative and therefore contributes with a trivial charge. The next 
orders in A yield the following equations 



d x f -(d x \nq + 2f 1 )f = 2f 1 , 
dxf-i ~ (d t hiq + 2/0 /_! = ft + 2/ - qr, 
d x f_ 2 - (d x \nq + 2/0 /_ 2 = 2/ /_i + 2/_ 1; 
9 x /_ 3 - (^lng + 2/0 /_ 3 = f_ x + 2/ /_ 2 + 2/_ 2 , 



(5.16) 
(5.17) 
(5.18) 
(5.19) 
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Integrating (5.16) we obtain 



(5.20) 



which again generates a trivial charge. From (5.17) we have 

2 N 



f- 



Q 



q 



idfq-d- 1 (r (d^qf 



(5.21) 



which is extremely nonlocal and the next densities becomes even more cumbersome. At 
this point we should mention that there is in fact three sets of nonlocal charges. The first 



set of charges comes from continuing this procedure, integrating (5.18) and so on. The 



second set of charges comes from choosing a trivial solution of the Bernoulli equation, 
fi = 0, and the solution f = qof (5.16). The third set of charges comes from the trivial 
choices fi = and fo = 0. All the subsequent densities are determined by the initial 
choices of f\ and f . Therefore, let fi = and fo = q, then we have 



fo = q, 
f-i = qd" 1 (q-r + 2), 
f. 2 = 2qd- 1 ((l + q)d- 1 (q-r + 2)), 

f-3 = qd-' (q {d-\q-r + 2)) 2 ) 

+ Aqd- 1 ((1 + q)d7 x 1 {{l + q)d-\q 



2))) 



(5.22) 
(5.23) 
(5.24) 

(5.25) 



Let fi = and fo = 0, then 

/_! = -qd-\ (5.26) 
/_ 2 = -2qd~\ (5.27) 

/_ 3 = -Aqd^r + qd~ l (q (^V) 2 ) , (5.28) 
The respective conserved charges are given by 

POO 

H- n = / f. n dx, n = 1,2,3,... (5.29) 



Let us emphasize that the charges are conserved for the whole hierarchy of equations, as 
can be explicitly checked for some of the individual models, either from the positive or 
negative flows. 



15 



6 Gauge and reciprocal transformations 



It was proved [I] that ( 4.25[ ) does not have real valued, smooth, soliton like solutions in the 



form u = u(x + ct). The models within the WKI hierarchy possess the loop-soliton kind 
of solutions, and in this section we aim to explain the algebraic origin of this peculiarity. 

6.1 Dressing method 

First of all, we will show that the dressing method [3], well known to construct soliton 
solutions of integrable hierarchies, is unable to solve the WKI hierarchy. 

The dressing method reconstructs the general operator U = E^ + A^[(j)] by gauging 
a vacuum solution Uq = E^\ corresponding to the trivial field configuration = 0. There 
must exist two dressing operators Q± effecting the following gauge transformations 

u = GiTOi 1 - (d x e±) ei 1 . (6.i) 



From (5.1) we see that there is a gauge freedom for the operator U = — (d x ^) 1 i r "~ 1 , i.e. 
\& i — y tyh where h is a constant group element. Consider the spectral problem for the 
vacuum (d x + Uq) \1/o — 0. Then, the transformations \l/o ^ ©+^0 an d ^ ©-^o will 
reconstruct the same operator U if they are gauge equivalent, i.e. O+^o = Q-^foh. In 



other words, the relations (6.1) are valid, if and only if, the dressing operators are related 



through a Riemann-Hilbert problem 

0: 1 + = \S? h%\ (6.2) 
The dressing operators are further factorized by a Gauss decomposition 

+ = e xm e x(1> ■ ■ ■ , 9_ = e yC V ( - 15 - - • , (6.3) 



where g Qv\ We emphasize that (6.2) is essential to construct explicit solutions 



and this relation is valid only upon the existence of both transformations 0±. 



Now consider the WKI operator replaced into (6.1) 



H 1 + qEl + tE\ = Q^Qg 1 - (d x G±) Q±\ (6.4) 

Taking this transformation with + we see that X® = is a solution of the projec- 
tion. The projection yields 

q E 1 + + rE 1 _ = -d x e x(1) e- xW , (6.5) 

which determines X^ 1 ) in terms of q and r. We can solve the higher grade projections 
recursively in terms of the previous terms to conclude that, in principle, the problem is 
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solvable. This means that 0+ is able to reconstruct U . However, this is not enough. The 
transformation (6.4) with 0_ still needs to work. Considering its projection we have 



tfi + qE \ + rE\ = e Y{0) H l e- Y(0) . (6.6) 
Assuming the group parametrization 

B = e Y(0) = e* fi V ff V £ ° (6.7) 

we obtain 

H 1 + qE\ + rE l _ = (l + 2 X ipe 2 *) H 1 - 2tbe 2<t, El + (2 X + 2^x 2 e 20 ) Si. (6.8) 

This relation implies if) = or x — 0, consequently g = or r = showing that 0_ 
is unable to reconstruct U and therefore, the dressing method can not be employed to 
this case. Nevertheless, we can explore a gauge relation between the WKI and AKNS 
hierarchies 118, 191. 



6.2 AKNS solutions 

Consider the AKNS spectral problem in space-time coordinates (y, s) 

(d y + F) $ = 0, (d s + G) $ = 0, (6.9) 

where F and G are the zero curvature Lax pair. We are omitting the index n labelling 
each flow, i.e. G = G(n) and s = s n . Denoting the AKNS fields by w and z, the Lax 
operator F is given by 

F = H 1 + w(y, s)E+ + z(y, s)E°_. (6.10) 
From the pure gauge form wE+ + zE°_ — d y vc = — d y e x< ~ 0) e~ x ^\ where e x<0> = Be uc and 



B is given by the Gauss decomposition (6.7), we obtain the following relations [23J 

d y <P = xdyipe 2 ^ d s <P = *ljd sX e 24, } (6.11 

%X ~ X 2 d y 



w = -dyipe 2 *, z = -dyX - X 2 dy^. (6.12) 



The fields in B are directly related to the famous tau functions, obtained by projecting 



(6.2) between the highest weight states of the algebra (for more details we refer the reader 
to pi ESI). The adjoint relations are (E j + ) ] = EZ\ (H j ) ] = H~ j and & = c. The 
fundamental states are {|/io), \to)} obeying E±\m) = 0, = for j > 0, i = 0, 1 

and also c\(m) = = 0, H°\no) = and H°\fj,i) = Let us define the state 

I/X2) = E^lfii), thus the projection of the LHS of (6.2) yields 



(to 


0: 1 0+ 


to) 


(to 


0: 1 ©+ 


to) 



(to 


0: 1 0+ 


to) 


(to 


0: 1 ©+ 


to) 



^ = 1 — 1^-1^ 1 — r» X = ~, — 1 — r- (6.13) 



(to 


0: 1 0+ 


to) 


(to 


©: 1 0+ 


to) 
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From the RHS of (6.2) we define the tau functions, classified in terms of the group element 
h, through 

Tij = (^l^o^o Vi)> hJ = 0,1,2. (6.14) 



Therefore, 

= — , ip= — , X = — • (6-15) 

TOO r ll r ll 

The group element h generates soliton solutions when it is written in terms of vertex 
operators in the form 



A' 



h = n exp 



(6.16) 



i=l 



The vertex operators should satisfy eigenvalue equations, defining the dispersion relations 
for the hierarchy 

[Tj, yU + s n V {n)] = fi(y, s n )Ti, (6.17) 

where Uq = H 1 and Vo(n) = H n are the vacuum Lax pair. Thus, upon enforcing the 
nilpotency property of the vertices between the states, i.e. (fik \ (Ti) 2 = 0, we have the 
explicit form given by 

N 

tm = (jik\ II (l + e A ri) 



1=1 



(6.18) 



s k i + (ri)*^ 1 + • • ■ + (r 1 r 2 ) M e /i+/2 + • ■ ■ + (r 1 r 2 r 3 ) w e /i+/2+/3 + ■ ■ ■ 



The relations (6.15) are very important since they relate a general multi-soliton solution 
specified by (6.14) and (6.18), and in the following sections we will write all the results in 
terms of 0, ip and X- 

For the AKNS hierarchy we have two vertex operators 



J = -0O 



(6.19) 



where Ki is a complex parameter. These vertices satisfy the following eigenvalue equations 



[Ti,yU + s n V (n)} = +2 («# + (k^ s n ) I\ = r]i(y, s n )r i5 (6.20a) 
[T'^yUo + s n V (n)] = -2 (my + (/q)™ s n ) Ti = ^(y, s n )Ti, (6.20b) 

where nfZ correspond to a given flow of the hierarchy and we have defined the dispersion 
relations rfi and £j. For instance, considering h = e Fl e r ' 2 we obtain the most simple 
nontrivial solution 



1 . ^1^2 m +f 2 £2 

Too = 1 + 7 ^ e ' r i2 = e 42 , 

(Kl - K 2 )^ 

m = l + ? M \ 2 e^ 2 , r 21 = e*. 
(«i - K 2 ) 2 



(6.21) 
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Note that the tau functions have the same form for every model withing the hierarchy, 



only the power (/««)" in the dispersion relation (6.20) changes for different flows 



6.3 Gauge transformation 



Through a gauge transformation $ i— >■ $ = g$ the operator (6.10) is transformed into 

F = g(H 1 + wE° + + zE° ) g~ l - d y gg-\ (6.22) 
Let us assume that g is a solution of 

wE% + zE°_ = g~ x d y g, (6.23) 



thus g = B from (6.7) and the transformation (6.22) yields 



F = gB x g- x = (1 + 2x# 20 ) H l + (2ip + 2 X ^ 2 e 24> ) E l + - 2 X e 2 ^E 1 _. 
This operator is not in the form of the WKI operator 

U = H 1 + q{x, t)E\ + r(x, t)E x _ 



(6.24) 



(6.25) 



due to the coefficient of H 1 . If we impose an equality between (6.25) and (6.24) we fall 



into the same kind of inconsistency encountered in (6.8), implying ip = or \ = and 



then q = or r = 0. The correct way to "normalize" the operator (6.24), absorbing the 
coefficient of H 1 , is through a reciprocal transformation. 



6.4 Reciprocal Backlund transformation 

Let us use the results from [19] which asserts that, the continuity equation 

d s a(y,s) + d y (3(y,s) = 0, 
is transformed to the reciprocally associated equation 

da(x,t) + d x 5(x,t) = 0, 
by the transformation (y, s) H- (x, t) such that 

dx = ady — fids, dt = ds, 

together with its inverse 



provided that 



dy = ^dx — Sdt, ds = dt, 



1 (5 
7 = -, 6 = . 



a 



a 



(6.26) 

(6.27) 

(6.28) 
(6.29) 
(6.30) 
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This theorem is straightforward to prove, noting that (6.26) and (6.27) are consequences 



of the compatibility of mixed second derivatives of (6.28) and (6.29), respectively, while 



(6.30) ensures the reciprocity between both coordinate systems. 



Therefore, using the transformation (6.28) with the operator (6.24) at hands, we get 

(6.31) 



the following spectral problem in the coordinates (x, t) 

(d x + U)$ = Q, (d t + V)$ = 

where 



U 



a 



V = G + -F. 

a 



From (6.24) it is evident that defining a through 

a(y,8) = l+2 X iJ>e 2 * 



the spectral problem (6.31) is mapped into the WKI problem with the relations 

2^ + Ix^e 14, 



q(y,s) 



a 



r(y,s) 



2 X e 



2d> 



a 



(6.32) 



(6.33) 



(6.34) 



We still need to express the space-time dependence in the variables (x,t), so from (6.28) 
we have the implicit relation 



x(y,t) = J a(y, t)dy + const 
T12T21 



y + 2 



(6.35) 



dy + const. 



where we have used (6.33) and (6.15) to express it in terms of the tau functions. Writing 



also (6.34) in terms of the tau functions, 



q{y,s) 



, Tl2 ( (Too) 2 + T12T21 

r n I ( r oo) 2 + 2ti 2 t 2 i 



r{y,s) 



T11T21 



( r oo) + 2t 12 t 2 i 



(6.36) 



Equations (6.36) and ( |6.35[ ) are the general solution of all the models within the WKI 
hierarchy in terms of the AKNS tau functions, which are known. For instance, replacing 



(6.21) and dispersion (6.20) with n = 3, corresponding to the solutions of model (4.10) 



the integral in (6.35) is elementary (we also set the integration constant to zero) and an 



explicit solution is obtained, whose graph is sketched in Fig. [T] The mathematical origin 
behind the loop-soliton solution is an implicit relation between the coordinates, arising 
from a reciprocal transformation. Multi-loop-solitons can be constructed similarly using 



the general form (6.18) 
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-5 -2.5 -15 -10 -5 

x(y,0) x(y,0) 

Figure 1: Parametric plots of q(y,t) versus x(y,t) at instant t = 0. Similar 

behaviour is obtained for r(y,t). The solutions keep their form unchanged during 
t evolution. 



6.5 From SPE to sine-Gordon 



This section is a particular case of a more general result presented in the next section, 
but because the hodograph transformation (y, s) — > (x, t) previously introduced [Hi U\ E] 
comes out so naturally from our definition (6.33), we will consider it first and it will also 
serve as a motivation for the forthcoming section. The model (4.25) is obtained from 
(4.22) when v — u, or equivalently r = q, thus imposing this restriction upon (6.34) we 



get 



X 



1 + V 2 ' 



(6.37) 



which substituted in the first relation of (6.11) yields 



-2<t> 



X 



(6.38) 



Therefore, (6.33), note also (6.30), can be written in terms of one of the fields only 

1 - 



a = 7 



-i 



= 1 - x l 

i + i) 2 l + x 2 



2e 2 * - 1. 



(6.39) 



Reminding from the trigonometric identity cos 6 = 2 cos 2 | — 1, the relation (6.39 ) strongly 
suggests the parametrization 

1 = cos#, e 



a = 7 



cos 



~X = tan 



e 



(6.40) 



The transformation (6.28) implies u x = a l u y , but q = u x so from (6.34) we therefore 
have 



u. 



-2 X e 2<t) = 2 tan f cos 2 f = sin 9. 



(6.41) 
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Squaring u x = a 1 u y we have 



cos 



1 (1 



2 --"« cos 2 #) 



ft 



7 2 = 1 + u 2 x , 



(6.42) 



which is the relation from where [8] started. Deriving both sides of (6.42) and plugging 
into the equations of motion ( 4.22[ ) (with v = u) we have the conservation law (6.27) and 



its reciprocal (6.26), namely 

d a = d x (2m 2 7 ) , d s a = d y (-2u 2 ) 



(6.43) 



Using (6.40) and (6.41) in the second equation of (6.43), we finally get the sine-Gordon 
model 



u = \6 a , 



9 ys = 4sin#. 



(6.44) 



6.6 From 2-SPE to Lund-Regge 



Let us remind from the identities q = u x , r = v x , a = - and d x = a 1 d v . Taking the 



product of both terms in (6.34) yields 



7=1 + u x v x 



a" 



1 - UyVy. 



Deriving (6.45a) and replacing into (4.22a) one gets 



,77* u xV x t 



Vn 



v. 



+ 4w + 4 



UVJJ X uvu x v. 



- 2- 



v. 



v. 



+ 2u 2 x v + 2u ( 7 2 - 1) 



while from (4.22b) one has 

u x v xt 



u x v n uvu x v x 
4^— _)_ 2 — 



V r 



v x 



+ 2u 2 x v + 2u ( 7 2 - 1) 



(6.45a) 
(6.45b) 



(6.46) 



(6.47) 



Simplifying these two last equations we obtain the respective conservation law (6.26) and 



its reciprocal (6.27), which read 



dtl = d x (2uvy) , 
d s a = d y (—2uv) . 

Deriving (6.45b) with respect to s and using ( |6.48b ) we therefore have 

v y (u ys — Aau) + u y (v ys — 4ftt>) = 0. 
A consistency condition, e.g. under the choice v = u, requires that 

4cra, 



Vy S = 4ft V. 



(6.48a) 
(6.48b) 



(6.49) 



(6.50a) 
(6.50b) 
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From (6.50) we have u and v in terms of Xi 4>i ^> once u y = aq and v y = ar are given from 

fields 

(6.51) 



(6.34). However, for our purposes, it is more convenient to introduce the auxiliary fields 

ijj = i/je rp , x = X^- 



Then, the relations (6.11) read 



i + xi> : 



and (6.34) yields 



v y = -2 X et 



Replacing (6.53) into (6.50) we obtain 



u 



d s ipe~ 



d s xe q 



2 (1 + x$) 



(6.52) 



(6.53a) 
(6.53b) 



(6.54a) 
(6.54b) 



and if we employ (6.51) with (6.15), the solution in terms of tau functions reads 

ToodgTu - T 12 d s T 00 



u(y, s) 
v(y, s) 



2r oTn 

r ll r 00^s r 21 — r 2ld s T 00 
2T00 r 00 + r 12 r 21 



(6.55a) 
(6.55b) 



Deriving (6.54) with respect to y and equating to (6.53), after using (6.52), we obtain 



the Lund-Regge model |26[ l2~7] 



^t-4x(l + ^)=0. 
1 + xip 



Xys 



(6.56a) 
(6.56b) 



Therefore, we have explicitly established an equivalence between the 2-SPE (4.22) and 
the Lund-Regge model (6.56). 



As a particular case, choosing x — we have 



Ipys + 



1-^ 



-Aip (l -^ 2 ) = 0, 



(6.57) 



which is the sine-Gordon equation (6.44) through the parametrization ip = sin 
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7 Conclusions 



We have proposed a higher grading construction for integrable hierarchies. When the 
algebra Q — A\ with homogeneous gradation is chosen, this construction yields the WKI 
hierarchy. The previous known models [6J constitute the positive flows. We have ex- 
tended the WKI hierarchy to negative flows, where the first negative flow yields the 
two-component field generalization of the short pulse equation [3]. Therefore, it is clear 
that in all these models underlie the same algebraic structure. Some novel integrable 
nonautonomous models were also proposed, mixing a positive with a negative flow. These 
integrable mixed models may have applications in nonlinear optics, specially concerning 
accelerated ultra-short optical pulses [22J. 

We have derived local and nonlocal charges from the Riccati form of the spectral 
problem. The nonlocal charges are usually more difficult to be obtained, and we have 
shown how they can be constructed through positive power series expansion in the spectral 
parameter. This method can be used to other integrable hierarchies as well. 

We have demonstrated that the dressing method is unable to solve the WKI hierarchy, 
which in fact, does not have the usual solitonic type of solution, travelling with constant 
velocity. Combining gauge and reciprocal Backlund transformations we have established 
the precise connection between the whole WKI and AKNS hierarchies, mapping each 
flow of both hierarchies. This gives a formal explanation of various relations scattered 
through the literature, for instance, the connection between the elastic beam and mKdV 
equations |28j . the relation between the Dym and KdV equations [29J and also the re- 
cent hodograph transformation [5] relating the short pulse equation to the sine-Gordon 
model. We considered explicitly this mapping between the first negative flows of the 
WKI and AKNS hierarchies, and demonstrated the correspondence between the novel 
two-component short pulse equation (4.22) and the Lund-Regge model (6.56), general- 
izing the previous hodograph transformation (1.3). In the same way that the notorious 
Miura transformation relates the KdV and mKdV hierarchies, which in fact is a gauge 
transformation, the combination of gauge followed by a reciprocal transformation yields 
more involved relations. We believe that our results show this explicitly, and the same 
approach can be applied to other integrable hierarchies following our construction. This 
relations will generate exotic solitonic solutions, like the loop-solitons, arising from an 
implicit space-time dependence. 

We have constructed solutions for the whole WKI hierarchy, writing them in terms of 
the known AKNS tau functions. We have motivated that the models within the positive 
flows of the WKI hierarchy describe larger amplitude solutions than the corresponding 
usual solitonic solutions, which is expected from the higher nonlinearity contained in 
the equations of motion. We believe that these models may have particular interesting 
applications that were not further explored. 
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Finally, we stress that our higher grading construction is general and other affine Lie 
algebras can be considered, which will give rise to novel integrable models. For instance, 
a n-component generalization can be considered with the algebra A n _i ~ s£ n . The 
technique of using gauge plus reciprocal transformations will be well suited to deal with 
these cases. 
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